Responses of leaf water potential and gas exchange 
to the precipitation manipulation in two shrubs on 
the Chinese Loess Plateau 


LI Yangyang!?", CHEN Jiacun’, AI Shaoshui', SHI Hui? 


! State Key Laboratory of Soil Erosion and Dryland Farming on the Loess Plateau, Institute of Soil and Water Conservation, 
Northwest A&F University, Yangling 712100, China; 

? Institute of Soil and Water Conservation, Chinese Academy of Sciences, Yangling 712100, China; 

3 School of Environmental and Municipal Engineering, Xi'an University of Architecture and Technology, Xi'an 710055, China 


Abstract: Regulation of leaf gas exchange plays an important role in the survival of trees and shrubs under 
future climate change. However, the responses of leaf water potential and gas exchange of shrubs in 
semi-arid areas to the precipitation alteration are not clear. Here, we conducted a manipulated experiment 
with three levels of precipitation, i.e., a control with ambient precipitation, 50% above ambient 
precipitation (irrigation treatment), and 50% below ambient precipitation (drought treatment), with two 
common shrubs, Salix psammophila C. Wang & C. Y. Yang (isohydric plant, maintaining a constant leaf 
water potential by stomatal regulation) and Caragana korshinskii Kom. (anisohydric plant, having more 
variable leaf water potential), on the Chinese Loess Plateau in 2014 and 2015. We measured the seasonal 
variations of predawn and midday leaf water potential (Ya and Wma), two parameters of gas exchange, i.e., 
light-saturated assimilation (4n) and stomatal conductance (g), and other foliar and canopy traits. The 
isohydric $. psammophila had a similar A, and a higher g than the anisohydric C. korshinskii under drought 
treatment in 2015, inconsistent with the view that photosynthetic capacity of anisohydric plants is higher 
than isohydric plants under severe drought. The two shrubs differently responded to precipitation 
manipulation. Ya, An and g were higher under irrigation treatment than control for $. psammophila, and 
these three variables and Yna were significantly higher under irrigation treatment and lower under drought 
treatment than control for C. korshinskii. Leaf water potential and gas exchange responded to manipulated 
precipitation more strongly for C. Rorshinskii than for $. psammophila. However, precipitation manipulation 
did not alter the sensitivity of leaf gas exchange to vapor-pressure deficit and soil moisture in these two 
shrubs. Acclimation to long-term changes in soil moisture in these two shrubs was primarily attributed to 
the changes in leaf or canopy structure rather than leaf gas exchange. These findings will be useful for 
modeling canopy water-carbon exchange and elucidating the adaptive strategies of these two shrubs to 
future changes in precipitation. 
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1 Introduction 


Precipitation is one of the most important factors controlling primary productivity in terrestrial 
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ecosystems (Huxman et al., 2004). Ecosystems in areas with the low mean annual precipitation 
are predicted to be most susceptible to projected changes of precipitation associated with climatic 
warming (de Dios et al., 2007). The wind-water erosion crisscross region on the Chinese Loess 
Plateau has the most intensive losses of soil and water in the world. Mean annual precipitation in 
this area is about 400 mm, but mean annual evaporation is >1500 mm (Tang et al., 1993), so 
water scarcity is the primary factor limiting plant growth and survival in this ecosystem. 
Large-scale coal and oil mining have exacerbated the shortage of water resources by depleting 
underground water sources, and some plants already have symptoms of dieback (Fan, 2007). 
Understanding the influences of projected variations of precipitation on the physiology, growth 
and productivity of local vegetation is therefore needed. 

Plants alleviate the threat of drought in a short term by regulating stomatal closure. Therefore, 
plants were classified into two drought strategies based on stomatal regulation i.e., isohydry and 
anisohydry (Tardieu and Simonneau, 1998). Isohydric plants maintain leaf water potential ( Year) 
by rapidly closing their stomata during drought, thus decreasing stomatal conductance (gs) and 
photosynthetic assimilation. In contrast, anisohydric plants decrease Year during drought, which 
may increase their vulnerability to hydraulic failure. Anisohydric plants often close their stomata 
later, thus maintaining higher gs and carbon assimilation under severe drought (McDowell et al., 
2008). This framework has been extensively used to explain the mortality of drought-induced 
trees (Adams et al., 2009; Hartmann et al., 2013; Sevanto et al., 2014). Some experimental results, 
however, contradict this elaborate framework. For example, daily gs decreased more for 
anisohydric sunflower than for isohydric maize under severe water deficit (Tardieu and 
Simonneau, 1998). gs and photosynthesis were similar at the end of a dry period for co-existing 
isohydric and anisohydric Mediterranean woody species (Quero et al., 2011), the same pattern has 
been reported in different Vitis cultivars (Lovisolo et al., 2010). These results contradict the 
traditional view that anisohydric species can maintain higher levels of gas exchange under severe 
drought. Martinez-Vilalta and Garcia-Forner (2016) found that isohydric Pinus edulis exhibited 
chronic embolism, whereas anisohydric Juniperus monosperma had very little embolism, refuting 
the hypothesis that anisohydric species are more vulnerable than isohydric species to hydraulic 
failure. Garcia-Forner et al. (2016) suggested that iso/anisohydry in the regulation of Year may be 
independent of the dynamics of leaf gas exchange or the degree of hydraulic or carbon limitations 
under drought, necessitating the measurements of leaf water potential and gas exchange together 
rather than only inferring gas exchange or hydraulic traits from plant iso/anisohydric behavior. 

The long-term acclimation of plants to limited water supply is important for their growth and 
survival. Experimental manipulation of precipitation has provided a way to perform controlled 
studies of acclimation under the long-term drought, thus improving our understanding of adaptive 
strategies for different plants under the projected conditions of future climate change (Weltzin et 
al., 2003; Huang et al., 2017). Acclimation to the long-term drought may involve the adjustment 
of plant morphology, physiology, anatomy, growth and/or carbon partitioning among organs. 
Photosynthesis is one of the most crucial physiological processes involved in acclimation (Chaves 
et al., 2002; Flexas et al., 2004), and the effects of the long-term drought on photosynthetic 
acclimation have been studied. For example, net rate of light-saturated assimilation (An), gs, 
mesophylic conductance of CO2, maximum rate of carboxylation and maximum rate of electron 
transport all decreased when predawn Year was reduced by the manipulation of throughfall in the 
evergreen Quercus ilex, but the functional relationships between these parameters and Year were 
not affected by water treatments (Limousin et al., 2010). In contrast, An and gs for Q. ilex from a 
mesic habitat decreased more rapidly in mesic than in dry habitats (Martin-StPaul et al., 2012). An 
under full irrigation, specific leaf hydraulic conductance and Wear at zero assimilation in a mixed 
forest in southwestern USA decreased as precipitation decreased for both pifion and juniper trees, 
and leaf gas exchange acclimated to the level of precipitation in both species (Limousin et al., 
2013). The long-term water stress can also invoke acclimatory responses that decrease the 
sensitivity of photosynthetic capacity to drought in xeric but not in riparian Eucalyptus species 
(Zhou et al., 2016). These variations in acclimatory responses suggest that functional 
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photosynthetic acclimation may be influenced by the climatic and soil conditions but can also 
vary among species/populations. 

Salix psammophila C. Wang & C. Y. Yang (Salicaceae) and Caragana korshinskii Kom. 
(Fabaceae) are two typical deciduous shrubs in the wind-water erosion crisscross region on the 
Chinese Loess Plateau. Both shrubs play important roles in controlling soil erosion and blocking 
the invasion of sand. S. psammophila is important to the energy and paper industries, whereas C. 
korshinskii is economically important as a highly productive forage shrub. Effects of the 
short-term drought on seedling water relations and photosynthesis on these two species have been 
studied (Liu et al., 2007; Xu et al., 2012), but little is known about the responses of mature plants 
to the long-term drought and varied precipitation. Thus, we conducted a manipulated experiment 
with three levels of precipitation (a control with ambient precipitation, 50% above ambient 
precipitation (irrigation treatment), and 50% below ambient precipitation (drought treatment)) on 
these two shrubs to address the following three questions: (1) are An and gs higher in C. 
korshinskii than in S. psammophila under severe drought due to the anisohydric behavior of C. 
korshinskii and isohydric behavior of S. psammophila (Li et al., 2016)? (2) does the 
shallow-rooted S. psammophila respond more strongly than C. korshinskii to the higher 
precipitation, and is it severely inhibited under drought treatment? and (3) xeric species 
photosynthetically acclimate to the long-term drought better than riparian species (Martin-StPaul 
et al., 2012; Zhou et al., 2016), so does the drought-resistant C. korshinskii (Ma et al., 2008; Xu et 
al., 2012) photosynthetically acclimate better than the drought-vulnerable S. psammophila? 


2 Materials and methods 


2.1 Study area 


The study area was conducted in the Liudaogou Watershed of the Shenmu County, Shaanxi 
Province, China (38°46'—-38°S1'N, 110°21'-110°23'E; 1081.0-1273.9 m a.s.l.), which is in a 
typical wind and water erosion crisscross region on the Chinese Loess Plateau. The mean annual 
precipitation is 437 mm, with 70% falling between July and September. The annual mean 
temperatures is 8.4°C, with an annual accumulative temperature >10°C of 3248°C, 153 frost-free 
days, total annual solar radiation of 5922 MJ/m? and mean annual sunshine duration of 2836 h 
(Tang et al., 1993). The area has mainly loamy, sandy and silt-dammed soils. Vegetation includes 
plants common in loessial and desert regions; typical woody plants are mainly trees, such as 
Populus simonii and Salix matsudana and shrubs, such as S. psammophila, C. korshinskii, 
Amorpha fruticosa and Hedysarum scoparium. 


2.2 Experimental design 


Two adjacent stands of healthy and uniformly growing S. psammophila and C. korshinskii on the 
flat top of a slope in the watershed were selected as experimental sites. Plants at the sites had been 
planted 30 years ago at the densities of 1600 and 2100 individuals/hm? for S. psammophila and C. 
korshinskii, respectively. The soil at the sites was classified as typical sandy soil. The average 
canopy size and height were 3.71 (+0.14) and 3.04 (+0.07) m for S. psammophila and 1.78 (+0.08) 
and 2.24 (40.21) m for C. korshinskii, respectively. 

The experimental design consisted of nine, randomly arranged 50 m° plots, with three 
replicates of three treatments per stand. The treatments were: (1) an ambient control that received 
a normal amount of precipitation, (2) a drought treatment with 50% less precipitation, and (3) an 
irrigation treatment with 50% more precipitation. We selected these lower and higher amounts of 
precipitation based on mean annual precipitation in the last 30 a and on the precipitation gradient 
(200-600 mm) within the distribution of the two shrubs on the Chinese Loess Plateau. For the 
drought treatment, 50% of the precipitation was diverted by transparent V-shaped troughs 
mounted on a steel frame at an interval of 20.0 cm. The steel frame was 2.5 m above the canopy 
and was supported by ten steel posts. Each post was inserted 0.5 m into the soil and fixed by 
concrete. The intercepted precipitation was collected and diverted by a PVC trough beside the 
frame and was stored in a large tank and filtered for use as the irrigation water. The 
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photosynthetic photon flux density at the canopy position was generally 90%—95% of that in the 
control, and ground-level temperature was 1°C—4°C lower during the growing season. For the 
irrigation treatment, water was delivered from the tank using a pump driven by a solar panel after 
a rainfall event with eight equally spaced sprinklers mounted about 2.5 m above the canopy. 
Precipitation <5 mm was not collected by the troughs, so only approximately 45% of the 
precipitation was diverted for application under the irrigation treatment. Vertical asbestos sheets 
were buried to a depth of 60 cm between plots to restrict the lateral movement of water in shallow 
soil. Three neutron-probe access tubes were installed to a depth of 3 m in each replicate for 
periodically monitoring soil moisture. The experimental setups were installed on 11 May 2014 for 
C. korshinskii and 21 May 2014 for S. psammophila, and no rain fell between the two 
installations. 

Three representative plants centrally located within each plot were selected as target shrubs for 
measuring physiological parameters. Year, gas exchange and leaf area index (LAI) were 
determined for these 27 plants of each species during the growing seasons in 2014 and 2015, and 
foliar dry mass per unit area and foliar 5'3C were measured only during the growing season in 
2015. 


2.3 Environmental variables 


Environmental variables were continuously monitored at a nearby field microclimatic station. 
Solar radiation (Rs) was measured using an Apogee20 PRY-P pyranometer (CS300, Apogee 
Instruments Inc., Logan, USA), air temperature and relative humidity (RH) were measured using 
a temperature and RH probe (HMPI55A, Vaisala Com., Helsinki, Finland), wind speed was 
measured using an anemometer (Met One, Met One Instrument, Inc., OR, USA) and precipitation 
was measured using a tipping-bucket rain gauge (TE525MM, Texas Electronics, Texas, USA). 
Data were recorded on CR1000 data loggers (Campbell Scientific, Shepshed, UK) every 10 s, and 
averages or totals were stored every 10 min. Vapor-pressure deficit (VPD, kPa) was calculated as 
the difference between saturated vapor pressure (es) and actual vapor pressure (ea), where es (Pa) 
was calculated using the Murray's formula (Murray, 1967): 


e;=610.7 x 107-5%237.3) (1) 
€a=RHxe,/100, (2) 
where fis the temperature (°C); and RH is the relative humidity (%). 


2.4 Measurements of eat and soil moisture 


Predawn (Ypa) and midday (Yma) Wear values were measured on sunny days at about 20 d 
intervals during the growing seasons in 2014 and 2015 using a pressure chamber (PMS 1000, 
PMS Instruments, Corvallis, USA). Two to three healthy mature leaves on current-year twigs 
from the south-facing side of each plant (6—9 leaves per plot) were measured on each measuring 
date. Leaves for the Ypa and Yma measurements were collected on the same days at 05:00—06:00 
(LST, local standard time) and 12:00—14:00, respectively. 

Soil moisture at a depth of 0.2-3.0 m was measured in each block at 20-cm intervals using a 
neutron probe (CNC 503, Beijing Chaosheng Technology Co. LTD, Beijing, China) 
approximately every three weeks. The neutron probe was calibrated in situ by saturating the 
ground near the experimental sites and periodically measuring soil moisture using both neutron 
probe and gravimetric method at different depths. The relationship between probe and gravimetric 
measurements was used to infer soil moisture. Soil surface (O—20 cm) moisture was measured by 
gravimetric method. 


2.5 Measurements of gas exchange, mass per unit area and 5'°C of leaves 
We measured two parameters of leaf gas exchange, An and gs, using a portable gas exchange 
system (LI-6400, Li-Cor Inc., Lincoln, USA) on two mature leaves of current-year twigs from the 


south-facing side of three plants per block and species and on the same dates that Wear was 
measured. Ay and gs were measured from 09:00 to 12:00 on a sunny day, with the light source set 
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to a saturating photosynthetic photon flux density of 1800 umol/(m?-s), and ambient temperature, 
relative humidity and ambient CO2 concentration were 7.4°C—27.2°C, 14.1%—62.8% and 376—402 
umol/mol, respectively, depending on the measurement date and time of the day. At the end of 
each measurement, the leaves enclosed in the chamber were collected for determining their 
projected areas that were imaged on a flatbed scanner and measured using ImageJ software (US 
National Institutes of Health, Bethesda, USA) 

After the measurement of gas exchange, samples of leaves were pooled on 29 May and 24 
August 2015 with the other mature leaves collected from current-year twigs for measuring leaf 
mass per unit area (LMA) and foliar 5'°C. The leaves were oven-dried at 75°C for 48 h for 
determining dry mass. LMA was calculated as the ratio of dry mass to projected area. Ground 
dried leaves were analyzed for 5'°C using a continuous-flow isotope ratio mass spectrometer 
(DeltaPlus, Finnigan, Germany). 5'3C was calculated as follows: 

5 13C=(Rsample/Rstandara—1 )x 1 000, (3) 
where R is the isotope ratio (13C/'7C); and Rsample and Rstandara are the isotope ratios of the sample 
and standard, respectively. The isotope values are expressed in delta notation relative to the 
VPDB (Vienna Peedee Belemnites) standard. The standard deviation for replicate analyses was 
+0.3%o. 


2.6 LAI 


LAI was monitored for three plants in each replicated block using a LAI-2200 Plant Canopy 
Analyzer (LI-COR Biosciences, Lincoln, USA). Three measurements from three directions, at 
intervals of 120°, were taken at half of projected canopy size between 05:30 and 08:00 on a sunny 
day. These measurements were taken three times in each growing season. 


2.7 Data analyses 


Effects of treatments on soil moisture, Ypa, Wma, An, gs, LMA, 5'°C and LAI on each sampling 
date were tested using one-way analysis of variance (ANOVA), with individual plots as a 
replicate (n=3) for each species. A Tukey post hoc test was used to identify differences among 
treatments when an ANOVA identified a significant effect. The differences between species of 
Poa, Mina, An and gs for each treatment during the growing season in 2015 were detected using an 
independent t-test. Analysis of covariance (ANCOVA) was used to compare the regressions 
between gs and An, between VPD and An or gs, and between Ya and An or gs among treatments for 
each species. We modeled the sensitivities of the gas exchange parameters to VPD and Ypa as 
follows: 

y=-—mlnVPD+n, (4) 
where m is the sensitivity of An or gs to VPD; and n is the corresponding value of Ay or gs at 
VPD=1 kPa (Oren et al., 1999), and 

y=ae™, (5) 

where x is the Ypa; b is the sensitivity of Anor gs to Ya; and a is the corresponding value of An or 
gs when x is close to 0 (Zhou et al., 2016). Standard major axis was used to compare the 
difference of m or b between species using SMATR (Standardised Major Axis Tests and Routines) 
2.0 software (Falster et al., 2006). 


3 Results 


3.1 Climatic variables and soil moisture 


Mean annual precipitation was 439.2 mm in 2014 and 388.4 mm in 2015, and precipitation during 
the growing season was 373.3 mm in 2014 and 305.6 mm in 2015 (Fig. 1). When the precipitation 
values during the growing season and for the entire year in 2014 and 2015 were compared with 
the corresponding values of the multi-years (1980-2014), we found that precipitation values were 
almost the same as means of multi-years in 2014, but 65 mm less for the growing season and 45 
mm less for the entire year in 2015, indicating that 2014 was a normal year but 2015 was a 


ChinaXiv@ (FEAF! 


JOURNAL OF ARID LAND 
50 — 4 35 
(a) 2014 — VPD Solar radiation} (b) 2015 
| == Precipitation | uh 30 g 
40 | T j 25 ef 
= NATL ah AL sk | | ij A h 25 E 
l N LAWN I Pa 
= 30 yi via | | 1 z 

& | | \ W f IIA ly | i Zr is X 
20 | S|. $ 
3 > E 
~ io 0 5 
N 


f Lt thi es 


l May 1Jun IJul 1 Aug 1Sep 1Qct 1 May lJun 1 Jul 1 Aug 1 Sep 1 Oct 
Date Date 


Fig. 1 Seasonal variations in precipitation, solar radiation and vapor pressure deficit (VPD) during the growing 
seasons in 2014 (a) and 2015 (b) 


relatively dry year. Mean daytime VPD was higher in 2015 (1.46 kPa) than in 2014 (1.24 kPa), 
but mean daily solar radiation (20.7 MJ/(m?.d) in 2014 and 20.3 MJ/(m?.d) in 2015) did not 
significantly differ between the two growing seasons. 


Soil moisture for all treatments was lower in 2015 than in 2014 for both species (Fig. 2), as 
expected. Soil moisture was relatively low under all treatments at the end of June and August in 
both years, corresponding to the low precipitation during these periods. The treatments did not 
significantly affect soil moisture at the end of June and August for both species. However, soil 
moisture significantly changed in the other months. Soil moisture for S. psammophila was 
generally lower under drought treatment and higher under irrigation treatment than the control 
during the growing season in 2014 and rainy season in 2015, no difference of soil moisture among 
treatments at the end of growing season in 2014 and dry season in 2015 may be related with leaf 
except in June-August in 2014 when soil moisture under drought treatment was not significantly 
decreased compared to the control. The result might be explained by the fact that the experimental 
setup was installed in the near term. 

Salix psammophila Caragana korshinskii 


194 (a) 2014 —o-— Control | (b) 2014 
—#- Irrigation 
—4— Drought 


Soil moisture (%) 


Soil moisture (%) 
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Fig. 2 Variations of soil moisture within 1.0 m soil depth for Salix psammophila (a and c) and Caragana 
korshinskii (b and d) under different precipitation treatments in 2014 and 2015. Bars are standard errors, n=3. * 
and ** indicate significance among three treatments at P<0.05 and P<0.01 levels, respectively. 
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3.2 Leaf water potential 


Paa for S. psammophila during the growing seasons in 2014 and 2015 generally ranged from —0.1 
to —0.5 MPa, except for at the end of June, when Ya was the lowest (—0.6 MPa in 2014 and —0.7 
MPa in 2015, Figs. 3a and b). Ya for S. psammophila was significantly affected by the treatments 
on two of the nine sampling dates in 2014 (P<0.05; Fig. 3a) and on five of the eleven sampling 
dates in 2015 (P<0.05; Fig. 3b). Y%a for S. psammophila usually responded more to irrigation than 
drought treatment, especially on the dates when Ya was significantly affected by the treatments 
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Fig. 3 Variations in predawn ( Y%a) and midday ( Yma) leaf water potentials for Salix psammophila (a, b, e and f) 
and Caragana korshinskii (c, d, g and h) under different precipitation treatments during the growing seasons in 
2014 and 2015. Bars are standard errors, n=3. * and ** indicate significance among three treatments at P<0.05 
and P<0.01 levels, respectively. 
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(Figs. 3a and b; Table 1). Ya for C. korshinskii ranged from —0.3 to —1.3 MPa during the two 
growing seasons, and Ypa for C. korshinskii under control and irrigation treatments was relatively 
stable throughout the growing seasons in both years, but Ypa under drought treatment 
significantly decreased during the late growing season in both years (Figs. 3c and d). Seven of the 
twelve Ypa measurements in 2014 and all ten of the measurements in 2015 were significantly 
influenced by the treatments (P<0.05). Ya was significantly lower under drought treatment and 
higher under irrigation treatment during the growing seasons in 2014 and 2015, especially on the 
dates when Ya was significantly affected by the treatments (Figs. 3c and d; Table 1). 

Wna for S. psammophila ranged from —0.6 to —1.6 MPa during the two growing seasons and 
decreased during the late growing season in 2015 (Figs. 3e and f). Three of the nine Yma 
measurements in 2014 and four of the eleven measurements in 2015 were significantly affected 
by the treatments (P<0.05). Yma was lower under drought treatment in 2014 and higher under 
irrigation treatment in 2015 compared with the control on the dates when the treatments 
significantly affected Yma (Figs. 3e and f; Table 1). Yma for C. korshinskii ranged from —1.3 to 
—2.8 MPa during the two growing seasons (Figs. 3g and h). Three of the twelve Yma 
measurements in 2014 were significantly affected by the treatments (P<0.05). Nine of the ten Ying 
measurements in 2015 were significantly affected by the treatments (P<0.05), with the highest 
Wna under irrigation treatment, and the lowest under drought treatment (Fig. 3h; Table 2). 


Table 1 Number of days with significant treatment effects for irrigation and drought treatments when compared 
with the control for Salix psammophila and Caragana korshinskii during the growing seasons in 2014 and 2015 


Salix psammophila Caragana korshinskii 
Variable Irrigation Drought Irrigation Drought 
2014 2015 2014 2015 N 2014 2015 2014 2015 M 
Pa 1 (+) 4 (+) 1 (-) 1 (-) 9, 11 3 (+) 9 (+) 7 (-) 10 (-) 12, 10 
Yaa 0 4 (+) 3(-) 2 (-)/2 (+) 9, 11 2 (+) 8 (+) 2 (-) 9 (-) 12, 10 
Ay 0 3 (+) 1 (-) 2 (-) 8,8 0 6 (+) 0 6 (-) 10, 8 
8gs 2 (+) 3 (+) 0 1 (-) 8,8 0 4 (+) 0 5 (-) 10, 8 


Note: Ya and Yna were predawn and midday leaf water potentials, respectively, A, and gs were light-saturated net assimilation and 
stomatal conductance, respectively. N was the measuring times in 2014 and 2015. The positive (+) signs indicate the increases and the 
negative (-) signs indicate the decreases of the number of days. 


Table 2 Averaged leaf water potentials and gas exchange parameters for Salix psammophila and Caragana 
korshinskii during the growing season in 2015 


Salix psammophila Caragana korshinskii 
Variable 
Control Irrigation Drought Control Irrigation Drought 
Ya (MPa) —0.30 —0.27 (10.00) —0.34 (-13.33) -0.65" -0.50 (23.08)* —0.91 (-40.00)* 
Yna (MPa) —1.09 —0.97 (11.01) —1.09 (0.00) -2.11° —1.93 (8.53)° -2.51 (-18.96)" 
An (umol/(m”+s)) 13.85 16.71 (20.65) 12.16 (—12.20) 16.24° 21.66 (33.37)° 11.33 (—30.23)"* 
gs(mmol/(m?-s)) 119.00 142.50 (19.80) 106.10 (—10.85) 116.00" 159.20 (37.25) 75.80 (34.66) 


Note: Yoa and Yna were predawn and midday leaf water potentials, respectively, and were measured for 11 times in Salix psammophila 
and 10 times in Caragana korshinskii. Ay and g; were light-saturated net assimilation and stomatal conductance, respectively and were 
measured 8 times for both species. The numbers within parentheses are the percentage increase relative to the control value, with 
positive value indicating increase and negative value indicating decrease. * and " indicate significance between two species at P<0.05 
and P>0.05 levels, respectively. 

Poa and Ying were generally lower for C. korshinskii than for S. psammophila in both years 

P P. p 

(Fig. 3; Table 2). There was no significant correlation between Y%a and Yra for S. psammophila 
(Fig. 4a), but Ypa was positively correlated with the water potential gradient (A% Pha Yma, 
P<0.05; Fig. 4b). Ypa was positively correlated with Yma for C. korshinskii (P<0.05; Fig. 4c), but 
Poa was not significantly correlated with the water potential gradient (Fig. 4d). 


3.3 Leaf gas exchange 


Parameters of leaf gas exchange for S. psammophila seasonally fluctuated, with An and gs much 
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Fig. 4 Relationships of predawn leaf water potential ( Ypa) with midday leaf water potential ( Yma) and water 
potential gradient (AY% Ypa- Yma) for Salix psammophila (a and b) and Caragana korshinskii (c and d) under 
different precipitation treatments. Bars are standard errors. Data from 2014 and 2015 were pooled together, in 
which n=9 in 2014 and 11 in 2015 for S. psammophila, and 12 in 2014 and 10 in 2015 for C. korshinskii. 


lower at the end of June and in mid-to-late August during the two growing seasons (Fig. 5). These 
decreases corresponded to the lower precipitation and soil moisture during these periods (Figs. 1 
and 2). Only one An and two gs of the eight measurements of leaf gas exchange for S. 
psammophila in 2014 were significantly affected by the treatments (P<0.05; Figs. 5a and b). Four 
An and three gs of the eight measurements in 2015 were significantly affected by the treatments 
(P<0.05; Figs. 5c and d). An and gs responded more to the irrigation treatment on the sampling 
dates when the treatment had significant effects (Figs.5a and d; Tables 1 and 2). An for C. 
korshinskii peaked in mid-to-late August during the two growing seasons. Treatments did not 
significantly affect An or gs in 2014 (Figs. 5e and f). All eight An and gs measurements were 
significantly affected by the treatments in 2015 (P<0.05), with the highest Anand gs values under 
irrigation treatment, followed by the control and the lowest under drought treatment (Figs. 5g and 
h; Tables 1 and 2). 

An and gs for both species responded to precipitation manipulation more in 2015 than in 2014 
(Fig. 5). An and gs during the growing season in 2015 were higher for C. korshinskii than for S. 
psammophila under irrigation treatment, but An was similar and gs was even lower than that of S. 
psammophila under drought treatment (Table 2), indicating that C. korshinskii did not have a 
higher photosynthetic capacity than S. psammophila under severe drought. The response pattern 
and magnitude of An and gs to precipitation manipulation differed between the two species. An and 
gs for S. psammophila clearly responded to irrigation but seldom to drought treatment. In contrast, 
An and gs for C. korshinskii responded strongly to both irrigation and drought treatments relative 
to the control (Tables 1 and 2). The response magnitude of both irrigation and drought treatments 
was larger for C. korshinskii than for S. psammophila (Table 2). Precipitation treatments did not 
affect the linear relationship between A, and g, for these two species (Fig. 6). 


3.4 Sensitivities of leaf gas exchange to VPD and Ypa 
Treatment and VPD or Ya did not significantly interact with An and gs for either species (P>0.05) 
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Fig.5 Variations in light-saturated net assimilation (An) and stomatal conductance (gs) for Salix psammophila (a, 
b, c and d) and Caragana korshinskii (e, f, g and h) under different precipitation treatments during the growing 
seasons in 2014 and 2015. Bars are standard errors, n=3. * and ** indicate significance among three treatments at 
P<0.05 and P<0.01 levels, respectively. 


in the ANCOVA, so the relationships between An and VPD or Ya, gs and VPD or Ypa did not 
differ among the treatments. When we pooled together the data from all treatments, leaf gas 
exchange decreased with increasing VPD and Ya, and decreased faster for S. psammophila than 
for C. korshinskii (P<0.05; Figs. 7 and 8), indicating a higher sensitivity of gas exchange to VPD 
and Yu. An and gs were the highest at Ya= —0.2 MPa for S. psammophila and at Ya= —0.4 MPa 
for C. korshinskii. Ya at a 50% loss of maximum An was —0.39 MPa for S. psammophila and 
—1.25 MPa for C. korshinskii, with a corresponding Ya at a 50% loss of maximum gs of —0.42 
MPa for S. psammophila and —1.21 MPa for C. korshinskii. Ppa corresponding to 50% losses of 
maximum An and gs was lower for C. korshinskii than for S. psammophila, suggesting that the 
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photosynthetic apparatus was more resistant to drought for C. korshinskii than for S. 
psammophila. 
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Fig. 6 Relationship between light-saturated net assimilation (An) and stomatal conductance (gs) for Salix 
psammophila (a) and Caragana korshinskii (b) under different precipitation treatments. Bars are standard errors. 
Data from 2014 and 2015 were pooled together, in which n=48 for S. psammophila and 54 for C. korshinskii. 
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Fig. 7 Relationships of light-saturated net assimilation (An) and stomatal conductance (gs) with VPD (vapor 
pressure deficit) for Salix psammophila (a and b) and Caragana korshinskii (c and d) under different precipitation 
treatments. Bars are standard errors. Data from 2014 and 2015 were pooled together, in which n=48 for S. 
psammophila and n=54 for C. korshinskii. 


3.5 LMA, 8C and canopy LAI 


Leaf 8°C and LMA for S. psammophila during the growing season in 2015 were significantly 
affected by the treatments (P<0.05), with 5!°C and LMA lower under irrigation treatment than the 
control. Treatments, however, did not significantly affect 6'°3C or LMA for C. korshinskii (Table 
3). Drought significantly decreased LAI for S. psammophila in 2015 and for C. korshinskii in both 
two seasons, and irrigation increased LAI for C. korshinskii during the middle growing season in 
2015 (Fig. 9). 
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Fig. 8 Relationships of light-saturated net assimilation (An) and stomatal conductance (gs) with predawn leaf 
water potential ( Ypa) for Salix psammophila (a and b) and Caragana korshinskii (c and d) under different 
precipitation treatments. Bars are standard errors. Data from 2014 and 2015 were pooled together, in which n=48 
for S. psammophila and n=54 for C. korshinskii. 


Table 3 Leaf mass per area (LMA) and 8'C for Salix psammophila and Caragana korshinskii during the 
growing season in 2015 


, LMA (g/m?) 33C (%0) 
Species Treatment 
May August May August 
Control 74.33+2.03* 123.05+5.40* —24.8+40.1* —25.6+0.2 
Salik Irrigation 61.1143.37° 99.1149.85° —25.4+0.1° ~26.040.1° 
psammophila 
Drought 72.00+2.08°* 125.09+410.7* —25.0+0.3 —25.2+0.1* 
Control 63.37+1.68 104.30+2.79 -—25.2+0.1 —26.2+0.4 
Caragana Irrigation 61.77+2.92 104.73+1.37 -25.5+0.4 -25.7+0.3 
korshinskii 
Drought 59.42+1.23 89.92+3.48 —24.9+0.2 —25.9+0.1 


Note: Mean+SE, n=3. Different lowercase letters indicate significance among different treatments at P<0.05 level. 


4 Discussion 


4.1 Strategies of stomatal regulation and carbon assimilation in two shrubs 


S. psammophila is an isohydric plant, whereas C. korshinskii is an anisohydric plant, based on 
seasonal changes in Year, daily stomatal regulatory capacity and stem hydraulic vulnerability 
(Dong and Zhang, 2001; Jiang and Zhu, 2001; Xu et al., 2012; Li et al., 2016). Seasonal 
fluctuations in Ya and Yna for S. psammophila were small, and Ya and Yma were not 
significantly affected by drought and correlated for most sampling dates (Fig. 3). Ypa and water 
potential gradient were significantly correlated, suggesting that S. psammophila could 
maintain %a and adjusted its water potential gradient when soil water potential decreased, 
behaving as a typical isohydric plant (Tardieu and Simonneau, 1998; Martinez-Vilalta and 
Garcia-Forner, 2016). In contrast, seasonal fluctuations in Ypa and Yna and the effects of drought 
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Fig. 9 Leaf area index of Salix psammophila (a) and Caragana korshinskii (b) under different precipitation 
treatments during the growing seasons in 2014 and 2015. Bars are standard errors, n=3. Different lowercase letters 
indicate significant differences among different treatments at P<0.05 level. 


on Ya and Yna were larger for C. korshinskii than for S. psammophila (Fig. 2). The linear 
regression between Ya and Yna was significant, but Ypa was not significantly correlated with the 
water potential gradient for C. korshinskii, suggesting that Yma decreased as soil water potential 
decreased and that this shrub behaved as an anisohydric plant (Tardieu and Simonneau, 1998; 
Martinez-Vilalta and Garcia-Forner, 2016). 

An and gs were higher for C. korshinskii than for S. psammophila under irrigation treatment, 
perhaps due to the higher leaf N content (1.84% for C. korshinskii and 1.05% for S. psammophila) 
of this leguminous species. An was similar for C. korshinskii and S. psammophila, and gs was 
lower for C. korshinskii than for S. psammophila under severe drought (Table 2), inconsistent 
with the higher An and gs for anisohydric than for isohydric plants under extreme drought. This 
pattern has been reported for various crops, woody species and Vitis cultivars (Tardieu and 
Simonneau, 1998; Lovisolo et al., 2010; Quero et al., 2011), further suggesting that 
iso/anisohydric behavior based on the regulation of Year may be difficult to associate with their 
gas exchange behavior (Martinze-Vilalta and Garcia-Forner, 2016). 


4.2 Responses of Y-ar and gas exchange to precipitation manipulation 


The parameters of Wear and gas exchange for both shrubs responded to the precipitation 
manipulation, but the patterns of response differed between the two species. The response of S. 
psammophila to irrigation and drought treatments was asymmetrical. Pha, An and gs for S. 
psammophila were higher under irrigation treatment than under the control but were similar 
between drought and control treatments. Ya, Pma, An and gs for C. korshinskii all responded 
symmetrically to irrigation and drought treatments. All these variables were higher under 
irrigation treatment and lower under drought treatment than under the control (Figs. 3 and 5; 
Tables 1 and 2). These different patterns of response between the two shrubs may be due to the 
hydraulic constraints, especially root distribution and the acclimation of water uptake by roots. 

S. psammophila has shallow roots, with roots not exceeding 1.5 m, but the roots spread 
horizontally by up to 3-fold of the canopy size. This root configuration restricts S. psammophila 
to water supplies in shallow soil (Liu et al., 2010), so Year and gas exchange for this species are 
very sensitive to irrigation. Stronger stomatal regulation in the isohydric S. psammophila may 
have limited responses of Wear, An and gs to drought on most of the measuring dates. C. 
korshinskii has a taproot system with dense lateral roots, and roots can reach deeper than 6.0 m 
for 12—15-year-old plants (Niu et al., 2003). A hydrogen isotope experiment in the Ulanbuh 
Desert found that C. korshinskii mainly acquired water from depths of 0-30 and 60-90 cm during 
late autumn (Zhu et al., 2010), suggesting that its root system took up water from both shallow 
and deep layers. This versatile strategy for acquiring water likely accounts for the symmetrical 
patterns of the responses of Wear, An and gs to irrigation and drought treatments. The asymmetrical 
pattern of response of S. psammophila to the water treatments in our study was similar to the 
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pattern in the isohydric P. edulis in southwestern USA, whereas the symmetrical patterns for C. 
korshinskii differed from that for anisohydric J. monosperma (Limousin et al., 2013). The 
photosynthetic capacity of J. monosperma was similar in irrigation and control plots and was 
significantly lower under drought treatment (Limousin et al., 2013). These differences in response 
to precipitation manipulation between two anisohydric species may be due to a greater reliance on 
deep soil water in J. monosperma than in C. korshinskii. 

C. korshinskii responded more strongly than S. psammophila to the precipitation manipulation, 
outperformed S. psammophila under irrigation treatment and was severely inhibited under 
drought treatment (Table 2). Leaf N content was higher for C. korshinskii than for S. 
psammophila due to its function of fixing N, which may help to increase photosynthetic capacity 
under irrigation in the nutrient-limited area. The high infiltration capacity of the sandy soils at the 
sites may also have allowed rainwater to quickly infiltrate to deep soil, which would be difficult 
to access by the shallow-rooted S. psammophila. Our previous data for the growth of these two 
shrubs (Ai et al., 2017) also demonstrated that seasonal increase in aboveground biomass was 
larger for C. korshinskii than for S. psammophila under irrigation treatment but smaller under 
drought treatment, supporting the validity of our data for gas exchange. The greater plasticity of 
Year and gas exchange for C. korshinskii may increase its fitness in highly variable environments, 
partly contributing to the wide distribution of this species on the Chinese Loess Plateau. The 
stricter stomatal regulation and more conservative response of gas exchange to drought for S. 
psammophila, however, may increase its adaptation in drought and desert areas, so the two shrubs 
may have different strategies for adapting to a long-term drought. 


4.3 Acclimation of leaf gas exchange to the long-term alteration of water supply 


The sensitivities of leaf An and g, to VPD or Ypa did not differ among the treatments for either 
species, suggesting that VPD and precipitation manipulation did not alter the functional 
relationship between leaf gas exchange and water stress. These findings are consistent with those 
reported for the evergreen trees Q. ilex and Arbutus unedo in a Mediterranean Macchia ecosystem 
(Nogués and Alegre, 2002; Ripullone et al., 2009; Limousin et al., 2010; Mission et al., 2010), 
suggesting that two years of variation in soil moisture did not acclimate photosynthetic gas 
exchange in the plants. Studies of Q. ilex across a precipitation gradient (Martin-StPaul et al., 
2012), pifion and juniper trees under precipitation manipulation in southwestern USA (Limousin 
et al., 2013) and xeric Eucalyptus species (Zhou et al., 2016), however, have reported acclimatory 
responses in the sensitivity of photosynthetic capacity to a long-term drought. These 
contradictions among studies may be due to the differences in species, developmental stage 
(seedling, sapling and mature tree) and degree of water stress. 

Long-term structural adjustment may alleviate the effects of drought on photosynthesis and gs 
(Maseda and Fernandez, 2006). The treatments of precipitation manipulation influenced LMA for 
S. psammophila, which was lower under irrigation treatment. The higher photosynthetic rate and 
lower integrated water-use efficiency (lower foliar 5'°C) for S. psammophila under irrigation 
treatment than under the control may therefore be partly attributed to the changes in foliar 
structural properties. Leaves with a lower LMA generally have higher rates of An and lower 
water-use efficiencies (Lamont et al., 2002; Gulias et al., 2003). LMA for C. korshinskii was not 
significantly influenced by the treatments, so the responses of leaf gas exchange to the irrigation 
and drought treatments were likely due to factors other than the alteration of foliar structure. LAI 
was significantly lower under drought treatment than under the control for both shrubs and was 
partly higher under irrigation treatment in 2015 (Fig. 8), indicating an acclimatory response of 
canopy structure. Changes to hydraulic conductance in soil-plant systems may be responsible for 
altered leaf gas exchange in precipitation manipulation (Ripullone et al., 2007; Limousin et al., 
2013; Zhou et al., 2016), so foliar and canopy properties may serve as potential regulators of 
water uptake by roots and hydraulic transport, possibly affecting photosynthesis in these two 
shrubs. 

Finally, this study identified the behavior of mature shrubs in a natural environment in response 
to projected variation in precipitation, but our experiments were conducted only for two years. 
Long-term experiments are thus needed to characterize the drought resistance and resilience of 
these fragile ecosystems to projected climate change. 


LI Yangyang et al.: Responses of leaf water potential and gas exchange to the precipitation manipulation... 


5 Conclusions 


Two consecutive years of precipitation manipulation did not alter the strategy of stomatal control 
in the two shrubs. An was similar for isohydric S. psammophila and anisohydric C. korshinskii, 
and gs was even higher for S. psammophila under severe drought, inconsistent with the view that 
An and gs under extreme drought are higher for anisohydric plants. Parameters of Wear and gas 
exchange ( Ypa, Pma, An and gs) for S. psammophila were most strongly influenced by irrigation 
treatment and for C. korshinskii were significantly influenced by both irrigation and drought 
treatments. Year and gas exchange were more plastic to irrigation and drought treatments for C. 
korshinskii. Precipitation manipulation did not alter the sensitivity of leaf gas exchange to VPD or 
soil moisture in either species. Acclimation to changes in soil moisture for the two species was 
primarily due to the changes in foliar structure (e.g., LMA) and/or canopy structure (e.g., LAT) 
rather than photosynthetic gas exchange. These findings may be useful for modeling water-carbon 
exchange in canopy and elucidating the adaptive strategies of these two shrubs to future changes 
in precipitation. 
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